Microscopic modelling of defects production 
and their annealing after irradiation in silicon 
for HEP particle detectors 



S. Lazanu^'\ I. Lazanu^ and M. Bruzzi*^ 

^National Institute for Materials Physics, POBox MG-7, Bucharest- Magurele, 

■ Romania 
O ■ 

' University of Bucharest, Department of Nuclear Physics, POBox MG-11, 

■ Bucharest- Magurele, Romania 

Q ' ^Universita di Firenze, Dipartamento di Energetica, Via S. Marta 3, 50139- 

C/^ ■ Florence, Raly 

(N 

> . 

CN . Abstract 



o 

0^ 



In this contribution, the production of defects in radiation fields and their evolution 



I toward equilibrium in silicon for detector uses has been modelled. In the quantita- 

CN ■ tive model developed, the generation rate of primary defects is calculated starting 

from the projectile - silicon interaction and from recoil energy redistribution in the 



lattice. Vacancy-interstitial annihilation, interstitial migration to sinks, divacancy 
and vacancy-impurity complex {VP, VO, V2O, CiOi and CiCg) formation are con- 
sidered. 

' The results of the model support the experimental available data. 

The correlation between the initial material parameters, temperature, irradiation 
^ ■ and annealing history is established. The model predictions could be a useful clue 

in obtaining harder materials for detectors at the new generation of accelerators or 
for space missions. 
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1 Introduction 



Silicon is the most used semiconductor material for detectors in particle physics 
experiments. The principal obstacles to long term operation in the extreme ra- 
diation conditions are the changes in detector parameters, consequence of the 
modifications in material properties during and after irradiation, ultimately 
related to the defects induced this way in the lattice. 

A point defect in a crystal is an entity that causes an interruption in the lattice 
periodicity. In this paper, the terminology and definitions in agreement with 
M. Lannoo and J. Bourgoin [1] are used in relation to defects. 

We denote the displacement defects, vacancies and interstitials, as primary 
point defects, prior to any further rearrangement. 

The microscopic model developed in this contribution is able to explain quanti- 
tatively, without free parameters, the production of primary defects in silicon, 
and their evolution toward equilibrium because in silicon, these defects are es- 
sentially unstable and interact via migration, recombination, and annihilation 
or produce other defects. 

The generation rate of primary defects is considered in a large range of values, 
and is calculated starting from the projectile - silicon nucleus interaction, and 
from the analytical extensions of the Lindhard theory [2]. 

Silicon used for detectors in high energy physics is n-type high resistivity 
(1 -i-6KQ ■ cm) phosphorus doped FZ material. The concentrations of intersti- 
tial oxygen Oi and substitutional carbon Cs in silicon are strongly dependent 
on the growth technique. In high purity Float Zone Si, oxygen interstitial 
concentrations are around 10^^ cm~^, while in Czochralski Si these concen- 
trations can reach values as high as 10^^ cm~^. Because Czochralski silicon is 
not available in detector grade quality, an oxygenation technique developed 
at BNL produces Diffusion Oxygenated Float Zone in silicon, obtaining a Oi 
concentration of the order 5 • 10^^ cm~^. These materials can be enriched in 
substitutional carbon up to [Cj] ^ 1.8 ■ 10^^ cm~^. 

Vacancy-interstitial annihilation, interstitial migration to sinks, divacancy and 
vacancy-impurity complex {VP, VO, V2O, CiOi and CiCs) formation are con- 
sidered. 

Model predictions are compared with experimental data. 

The correlation between the initial material parameters, temperature, irradi- 
ation and annealing history is established. 
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2 Production of primary defects and their kinetics 



The basic assumption of the present model is that the primary defects, va- 
cancies and interstitials, are produced in equal quantities and are uniformly 
distributed in the material bulk. They are produced by the incoming particle, 
or as a consequence of the subsequent collisions of the primary recoil in the 
lattice. 

Prior to the irradiation process, in silicon there are thermally generated defects 
(only Frenkel pairs are considered). 

After the irradiation, the following stable defects have been identified in silicon 
(see References [1,3]): Sii, VP, VO, V2, V2O, QOi, d and QC,. 

The pre-existing thermal defects and those produced by irradiation, as well as 
the impurities, are assumed to be randomly distributed in the solid. An im- 
portant part of the vacancies and interstitials annihilate. The sample contains 
certain concentrations of impurities, which can trap interstitials and vacancies 
respectively, and form stable defects. 

Vacancy-interstitial annihilation, interstitial migration to sinks, divacancy, va- 
cancy and interstitial impurity complex formation are considered. The role of 
phosphorus, oxygen and carbon is taken into account, and the following stable 
defects: VP, VO, V2, V2O, CiOi, Ci and CiCs are considered. Other possible 
defects as V3O, V2O2, V3O3 [4], are not included in the present model. 

The following picture describes in terms of chemical reactions the mechanisms 
of production and evolution of the defects considered in the present paper: 

V + I ^annihilation (1) 

G 



I ^ sinks (2) 



V + oSiVO (3) 

V + pS^VP (4) 

V + V^V2 (5) 

K3 

V + VO^ V2O (6) 

I + Cs'^a (7) 



VO is the A centre. 
VP is the E centre. 
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Q + CO, (8) 
A + J^O (9) 

Q + Cs"^ CiCs (10) 

The reaction constants f^j (i = 1, 3 -r 8) have the general form: 

K, = Cue^v{-E,lkBT) (11) 

with V the vibration frequency of the lattice, Ei the associated activation 
energy and C a numerical constant that accounts for the symmetry of the 
defect in the lattice. 

The reaction constant related to the migration of interstitials to sinks could 
be expressed as: 

K2 = av\^ exp {-E^lkBT) (12) 
with a the sink concentration and A the jump distance. 

The system of coupled differential equations corresponding to the reaction 
scheme (1)-^ (10) cannot be solved analytically and a numerical procedure 
was used. 

The following values of the parameters have been used: E\=E2 = 0.4 eV, E-^ 
= 0.8 eV, = 1.4 eV, E5 = 1.1 eV, Eq = 1.3 eV, E7 = 1.6 eV, Eg = 1.7 eV, 
u = 10^3 Hz, A = 10^5 cm2, a = 10^° cm'^. 

The generation term is the sum of two components: 

G = Gr + Gt (13) 

where Gr accounts for the generation by irradiation, and Gt for thermal 
generation. 

The concentration of the primary radiation induced defects per unit fluence 
(GPD) in silicon has been calculated using the explicit formula (see details, 
e.g. in references [5,6]): 

GPD {E) = ^ [y(^] L{En^sidn = -^^^^NIEL (E) (14) 

where E is the kinetic energy of the incident particle, Nsi is the atomic density 
in silicon, Asi is the silicon atomic number, Esi - the average threshold energy 
for displacements in the semiconductor, Ejn - the recoil energy of the residual 
nucleus produced in interaction i, L{ERi) - the Lindhard factor that describes 
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the partition of the recoil energy between ionisation and displacements and 
{da/dQ)i - the differential cross section of the interaction between the incident 
particle and the nucleus of the lattice for the process or mechanism i, respon- 
sible in defect production. Na is Avogadro's number. The formula gives also 
the relation with the non ionising energy loss {NIEL). 

Due to the important weight of annealing processes, as well as to their very 
short time scale, CPD is not a measurable physical quantity. 

If the anisotropy of the silicon lattice is not considered, the simplifying hypoth- 
esis of random distribution of CPD for all particles could be introduced: the 
identity of the particle is lost after the primary interaction, and two different 
particles could produce the same generation rate {Gr) of vacancy-interstitial 
pairs if the following condition: 

is fulfilled. Here, $ is the flux of particles (a) and (6) respectively, and Ei and 
E2 their corresponding kinetic energies. 

3 Results and discussion 

3. 1 Comparison between model predictions and experimental data 

The model predictions have been compared with experimental measurements. 
A difficulty in this comparison is the insufficient information in published 
papers regarding the characterisation of silicon, and the irradiation parameters 
and conditions for most of the data. 

It was underlined in the literature [7] that the ratio of VO to VP centres 
in electron irradiated silicon is proportional to the ratio between the concen- 
trations of oxygen and phosphorus in the sample. For electron irradiation, in 
Reference [8] a linear dependence of the V2 versus VO centre concentration has 
been found experimentally. In the present paper, the ratios of concentrations 
of V2 to VO centres and VO to VP ones has been calculated respectively in 
the frame of the model, for the material with the characteristics specified in 
Reference [8], and irradiated with 12 MeV electrons, in the conditions of the 
above mentioned article. The time dependence of these two ratios is repre- 
sented in Figures la and lb; annealing is considered both during and after 
irradiation. It could be seen that for the ratio of V2 and VO concentrations, 
the curves corresponding to different irradiation fluences are parallel, while 
the ratio of VO to VP concentrations is fluence independent, in the interval 
2 ■ 10^^ 5.5 • 10^"^ cm~^, in good agreement with the experimental evidence. 
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The ratio between V2 and VO concentrations is determined by the generation 
of primary defects by irradiation, while the ratio between VO and VP con- 
centrations is determined by the concentrations of oxygen and phosphorus in 
sihcon. 




Fig. 1. Time dependence for a)V2/VO and h)VO/VP concentrations calculated for 
silicon with 1.4 • 10^^ P/cm^, 5 • 10^^ O/cm^, 3 ■ 10^^ C/cm^, after 12 MeV electron 
irradiation, with the flux 5.8- 10^*^ e/cm^-s, up to the fluences: 2-10^^ e/cm^, 4-10^^ 
e/cm^, 8 • 10^2 e/cm^, 2 • 10^'' e/cm^, 2 • 10^^ e/cm^, 5.5 • lO^'' e/cm^ followed by 
relaxation (see reference [8]. 

Our estimations are also in agreement with the measurements presented in 
reference [9], after electron irradiation, where defect concentrations are pre- 
sented as a function of the time after irradiation. In Figure 2, both measured 
and calculated dependencies of the VO and VP concentrations are given. The 
irradiation was performed with 2.5 MeV electrons, up to a fiuence of 3 ■ 10^^ 
cm~^. The dependencies put in evidence the important role played by the 
carbon-related defects. The relative values are imposed by the arbitrary units 
of experimental data. 

A good agreement has also been obtained between the absolute values of 
concentrations of VP + V2 and CiCg centres predicted by the model, and 
the experimental results after 1 MeV neutron irradiation at a total fiuence of 
5.67 • 10^'^ cm~^, reported in reference [10]. The calculated 1.5 • 10^^ cm~^ and 
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Fig. 2. Time dependence of VO and VP concentrations after electron irradiation: 
crosses - experimental data from reference [9]; continuous line present model cal- 
culations; dashed line - without the consideration of carbon contribution to defect 
formation. 

4.1-10^^ cm~^ concentrations for VP + V2 and CiCg respectively, are in accord 
with the values of 1.1 • 10^^ cm~^ and 3.8 • 10^^ cm~^ , measured experimentally. 
For the VO concentration, a poorer concordance has been obtained. 

The formation and time evolution of stable defects depends on various factors, 
e.g. the concentrations of impurities pre-existent in the sample, the rate of 
generation, and the temperature and the irradiation history. 

3.2 The variation of the concentration of stable defects in respect to the ir- 
radiation rate and the influence of oxygen 

The effect of oxygen in irradiated silicon has been a subject of intensive stud- 
ies in remote past. In the last decade a lot of studies have been performed to 
investigate the influence of different impurities, especially oxygen and carbon, 
as possible ways to enhance the radiation hardness of silicon for detectors in 
the future generation of experiments in high energy physics - see, e.g. refer- 
ences [11,12]. Some authors consider that these impurities added to the silicon 
bulk modify the formation of electrically active defects, thus controlling the 
macroscopic device parameters. Empirically, it is considered that if the silicon 
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is enriched in oxygen, the capture of radiation-generated vacancies is favoured 
by the production of the pseudo-acceptor complex vacancy-oxygen. Interstitial 
oxygen acts as a sink for vacancies, thus reducing the probability of formation 
of the divacancy related complexes, associated with deeper levels inside the 
gap. 

Two types of silicon have been considered in the model calculations: "stan- 
dard" material, with the following doping concentrations: 10^^ cm~^ phospho- 
rous content, 2-10^^ cm~^ oxygen, and 3-10^^ cm~^ carbon, and "oxygenated" 
material, with the same concentrations of phosphorous and carbon, but with 
200 higher oxygen content, 4 • 10^'' cm~'^ . 

These two types of materials are considered to be exposed to continuous irra- 
diation, with Gr varying from 5 • lO^l^Jpairs/s, to 50 V^/pairs/s, from half to 
half order of magnitude, at 20°C, during 10 years. 

Thermal generation has been considered each time in the calculations. The 
results are illustrated in Figure 3. The highest generation rate corresponds 
here to a value 10 times higher than that estimated for the forward position 
in the tracker silicon detector at the new LHC hadron accelerator (determined 
from integration over the energy range of interest of the convoluted function 
between simulated hadron flux spectra [13] and the CPD), while the smallest 
generation rate corresponds to a value 10 times smaller than the rate produced 
by protons from cosmic rays, in an orbit near the Earth, at about 380 Km, as 
measured by the AMS Collaboration [14]. 

It could be seen that the presence of oxygen slows down the formation of all 
defects related to oxygen and vacancy, i.e. VO, VP, V2, V2O and CjOj, and 
increases the rate of formation of CjC^ centres. It could also be observed that 
the highest concentrations attained by the VO centres during this period, 
as well as those corresponding to the V2O ones increase, from "standard" to 
"oxygenated" silicon, in the same ratio as the corresponding oxygen content. 
Because the process of formation of both divacancies and V2O centres is, in 
an important measure, slowed down in oxygenated silicon, and the maximum 
values of divacancy concentrations are reduced with the increase of oxygen 
concentration, the generation current in the depleted regions of p-n junction 
detectors could be decreased by oxygen addition, especially at smaller rates 
than those associated with curves (1). 

3. 3 The effect of the temperature 

Both "standard" and "oxygenated" silicon are considered to be exposed to 
similar continuous irradiations, characterised by a generation rate Gn = 5 ■ 
lO^y/pairs/s, during ten years, at temperatures between -20°C and 20°C. 
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Fig. 3. Influence of irradiation rate on the time evolution of defect concentrations. 
The values of Gr are decreasing from 5 • 10^ VI pairs/s - curve (1) to 50 VI pairs/s 
- curve (17), from half to half order of magnitude. 

Thermal generation has been considered for each temperature. The results 
are represented in Figure 4 for the concentrations of VO, V2 and V2O centres 
respectively. While in " standard" silicon the decrease of the temperature slows 
down the rate of VO formation, a much weaker dependence has been found 
in silicon enriched in oxygen. 

For the V2O centres, the maximum concentrations correspond to the highest 
temperature, 20°C, and similar values are found at the end of the time in- 
terval considered, although the introduction rate of this defect has different 
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standard Si oxygenated Si 




Fig. 4. Time dependence of the concentrations VO, V2 and V2O centres, in the 
conditions of continuous irradiation with Gr = 5- 10'' VI pairs/s, in "standard" and 
"oxygenated" sihcon, for 20°C (red), 10°C (green), 0°C (blue), - 10°C (cyan) and 
-20°C (magenta) respectively. 



behaviours in the two types of silicon considered. 



In contrast to oxygen related centres, the V2 concentration is more sensitive to 
temperature in "oxygenated" silicon, where its values are much lower in respect 
to the corresponding ones for "standard" material. More, the concentration of 
V2 centres increases with the decrease of the temperature from 20 to -20"C in 
"oxygenated" silicon, and decreases in the "standard" material. 
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3-4 The correlation between the history of the irradiation process and defect 
concentrations 

In some previous papers [15,16], we demonstrated in concrete cases the im- 
portance of the sequence of irradiation processes, considering that the same 
total fluence can be attained in different situations: the ideal case of instan- 
taneous irradiation, irradiation in a single pulse followed by relaxation, and 
respectively continuos irradiation process. As expected, after instantaneous 
irradiation the concentrations of defects are higher in respect with "gradual" 
irradiation. 

After this analysis, the specific importance of the irradiation and annealing 
history (initial material parameters, type of irradiation particles, energetic 
source spectra, flux, irradiation temperature, measurement temperature, tem- 
perature and time between irradiation and measurement) on defect evolution 
must be underlined. 

Now, we present comparatively in Figure 5 the time evolution of stable defects 
in standard silicon, irradiated during 10'' s, when 5 • 10''^VI pairs/s have been 
generated (continuous line), and during 2-10'' s, with 2.5-10'^^/ pairs/s (dashed 
line), respectively. The values have been chosen to produce the same amount 
of vacancy - interstitial pairs. The behaviour of stable defect concentrations 
is analysed after irradiation. It could be seen that a shorter and stronger 
irradiation conduces to higher concentrations of V2 and VP centres in the 
asymptotic limit, with important consequences for the leakage current, due to 
the position in the band gap of silicon of these defect levels. 

Another scenario that could be of interest is the comparison of the effects 
of continuous and pulsed irradiation, with the same integral rate of primary 
defect production. The time evolution of the concentrations of stable defects 
is represented in Figure 6. The effects of continuous irradiation of standard 
silicon, at 20°C, with a rate of 2.5 ■ 10"^ VI pairs/s (continuous line) are com- 
pared with pulsed one: 10^ s. irradiation, with 5 • lO^VI pairs/s, and 10'^ s. 
relaxation. As in the previous case (Figure 5), the concentration of VP centres 
follows the best the irradiation type, while for the other defects the succession 
irradiation - relaxation is almost not observable. 



4 Summary and conclusions 

The production of defects in radiation fields, in silicon for detector uses, and 
their evolution toward equilibrium, has been microscopically modelled in the 
frame of a quantitative model, without free parameters. Vacancy-interstitial 
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Fig. 5. Time evolution of the concentrations of (up to down): VO, VP, V2, V2O, 
CiOi, and CiCs in "standard" Si, at 20°C for two types of irradiation: 10'' s. irra- 
diation with Gr = 5 • 10'' VI pairs/s (continuous hne), and 2 • 10'' s. Gr = 2.5 • lO'^ 
VI/s (dashed hne) respectively. 

annihilation, interstitial migration to sinks, divacancy and vacancy-impurity 
complexes (VP; VO; V2O, CiOi and CiCg) formation is considered. 

The model supports the experimental data and confirms the different studies 
performed to investigate the influence of oxygen in the enhancement of the 
radiation hardness of silicon for detectors. The VO defects in oxygen enriched 
silicon are favoured in respect to the other stable defects, so, for detector 
applications it is expected that the leakage current decreases after irradiation. 
At high oxygen concentrations, this defect saturates starting from low fiuences 
at high generation rates of defects. 

The generation rate of primary defects, calculated starting from the projec- 
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Fig. 6. Time evolution of the concentrations of (up to down): VO, VP, V2, V2O, 
CiOi, and CiCg for two types of irradiation: continuous line - continuous irradiation 
with 2.5 • 10''' VI pairs/s during 10^ sec, (continuous line), and rectangular pulses 
10^ s. irradiation with Gr = 5 • 10^ VI/s, and 10^ s relaxation. 

tile - silicon nucleus interaction and from the extended Lindhard theory, is 
considered in a large range of values. 

The obtained results suggest the importance of the conditions of irradiation, 
temperature and annealing history on defect kinetics. 
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